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Let {X, Y) be a random vector whose conditional excess probability 0{x, y) :— P{Y <y \ X > x) 
is of interest. Estimating this kind of probability is a delicate problem as soon as x tends to be 
large, since the conditioning event becomes an extreme set. Assume that {X, Y) is elliptically 
distributed, with a rapidly varying radial component. In this paper, three statistical procedures 
are proposed to estimate 9{x,y) for fixed x,y, with x large. They respectively make use of an 
approximation result of Abdous et al. (cf. Canad. J. Statist. 33 (2005) 317-334, Theorem 1), a 
new second order refinement of Abdous et al.'s Theorem 1, and a non-approximating method. 
The estimation of the conditional quantile function 6{x,-)^ for large fixed x is also addressed 
and these methods are compared via simulations. An illustration in the financial context is also 
given. 

Keywords: asymptotic independence; conditional excess probability; elliptic law; financial 
contagion; rapidly varying tails 

1. Introduction 

Consider two positively dependent market returns, X and Y . It is of practical importance 
to assess the possible contagion between X and Y . Contagion formalizes the fact that 
for large values x, the probability P{Y > y \ X > x) is greater than P{Y > y): see, for 
example, Abdous et al. [1] or Bradley and Taqqu [9, 10], among others. Besides, this 
conditional probability is also related to the tail dependence coefficient which has been 
widely investigated in the financial risk management context: see, for instance, Frahm 
et al. [22]. Therefore, the behavior of the conditional excess probability 0{x,y) := P{Y < 
y I AT > a;) is of practical interest, especially for large values of x. Estimating this kind of 
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probability is a delicate problem as soon as x tends to be large, since the conditioning 
event becomes an extreme set. "Large" here essentially means that x is beyond the largest 
value of the X observations so that the conditional empirical distribution function then 
fails to be of any use, even if the probability 9{x,y) in itself is not a small probability, 
nor close to 1. Alternative methods have to be considered. 

A classical approach is to call on multivariate extreme value theory. Many refined 
inference procedures have been developed, making use of the structure of multivariate 
max-stable distributions introduced by De Haan and Resnick [15], Pickands [37] and De 
Haan [14]. These procedures are successful in the rather general situation where {X,Y) 
arc asymptotically dependent (for the maxima) , which means hcuristically that X and Y 
can be simultaneously large (see, e.g., Resnick [38] or Bcirlant et al. [5] for more details). 

Efforts have recently been made to address the problem in the opposite case of asymp- 
totic independence. In some papers, attempts are made to provide models for joint tails 
(see, e.g., Ledford and Tawn [33, 34], Draisma et al. [18], Resnick [39], Maulik and Resnick 
[36]). In a parallel way, Heffernan and Tawn [28] explored modeling for multivariate tail 
regions which are not necessarily joint tails, and Heffernan and Resnick [27] provided a 
complementary mathematical framework in the bivariatc case. 

In these papers, the key assumption is that there exists a limit for the conditional 
distribution of Y , suitably centered and rcnormalizcd, given that X tends to infinity. 
This assumption was first checked for bivariate spherical distributions by Eddy and 
Gale [19] and Bcrman [8], Theorem 12.4.1. Abdous et al. [1] obtained it for bivariate 
elliptical distributions, Hashorva [25] for multivariate elliptical distributions, Balkcma 
and Embrechts [3] for generalized multivariate elliptical distributions and Hashorva et 
al. [26] for Dirichlet multivariate distributions. 

Elliptical distributions form a large family of multivariate laws which have received 
considerable attention, especially in the financial risk context; sec Artzner et al. [2], 
Embrechts et al. [20], Hult and Lindskog [30], among others. Assume from now on that 
{X, Y) is clliptically distributed; Theorem 1 of Abdous et al. [1] exhibits the asymptotic 
behavior of 9{x, y) when x ^ oo for such an elliptical pair. The appropriate rate y = y{x) 
is made explicit to get a non-degenerate behavior of ]in\x^ooO{x,y). This rate depends 
on the tail behavior of the radial random variable R defined by the relation = {X'^ — 
2pXY + Y'^)/{1 — p^), where p is the Pearson correlation cocfhcient between X and Y . 
The only parameters involved in y{x) and in the limiting distribution are the Pearson 
correlation coefficient p and the index of regular variation of R (say a) or an auxiliary 
function of R (say depending on whether R has a regularly or rapidly varying upper 
tail. 

In financial applications, the regularly varying behavior of the upper tails is commonly 
encountered. Abdous et al. [1] provided a simulation study in the specific case where 
R has a regularly varying tail. Existing estimators of p and a were used therein to 
obtain a practical way to estimate excess probabilities. However, this assumption fails to 
hold in some situations, as shown by Levy and Duchin [35] who compared the fit of 11 
distributions on a wide range of stock returns and investment horizons. They concluded 
that the logistic distribution, which has a rapidly varying tail, gives the best fit in most 
of the cases for weekly and monthly returns. 
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The aim of this paper is to focus on the case where the radial component R asso- 
ciated with the eUiptical pair {X,Y) has a rapidly varying upper tail. A second order 
approximation result is obtained, which refines Theorem 1 of Abdous et al. in the case 
of rapid variation of R. As we mentioned earlier, under the rapidly varying upper tail 
assumption, the obtained approximation involves an auxiliary function -0 which has to 
be estimated. To the best of our knowledge, estimation of the auxiliary function has 
never been considered in the literature. We propose three statistical procedures for the 
estimation of 9{x^ y) for fixed a;, y with x large. They respectively make use of Abdous et 
al.^s approximation result, its second order refinement and an alternative method which 
does not rely on an analytic approximation result. 

The conditional quantilcs - also called regression quantiles - of a response Y given a 
covariate X have received significant attention; see, for example, Koenker and Bassett 
[32] or Koenker [31]. More specifically, extremal conditional quantiles have proven to be 
useful in economics and financial applications (sec, e.g., Chcrnozhukov [11]). In the same 
spirit, one might be interested in estimating the conditional quantile function 9{x,-)'^ 
when X is fixed and large. This problem is also addressed in detail in the present paper 
and the simulation study performed is presented in terms of estimation of both 9 and 



The paper is organized as follows. The second order approximation result is presented 
in Section 2, as well as some remarks and examples illustrating the theorem. The sta- 
tistical procedures are described in Section 3, where a semi-parametric estimator of tp is 
proposed. Section 4 deals with a comparative simulation study, while Section 5 provides 
an application to the financial context, revisiting data studied by Levy and Duchin [35]. 
Some concluding comments arc given in Section 6. Proofs are deferred to the Appendix. 

2. Asymptotic approximation 

Consider a bivariate elliptical random vector {X,Y). General background on elliptical 
distributions can be found in, for example. Fang et al. [21]. One can focus, without loss of 
theoretical generality, on the standard case where EX = KY = and Var X = Var Y = 1. 
A convenient representation is then the following (see, e.g., Hult and Lindskog [30]): 
{X, Y) has a standard elliptical distribution with radial positive distribution function H 
and Pearson correlation coefficient p if it can be expressed as 



where R and U are independent, R has distribution function H with Ei?^ = 2 and U is 
uniformly distributed on [— 7t/2, 37t/2]. Hereafter, to avoid trivialities, we assume \p\ < 1. 
Let $ denote the normal distribution function and ip its density, that is. 



9{x,-r. 



{X,Y)^ i?(cos U, pcosU+ \/~V~f? sin U) , 




e 



This paper deals with elliptical distributions with rapidly varying marginal upper tails 
or, equivalently, with a rapidly varying radial component. More precisely, the radial 
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component R associated with (X, Y) is assumed to be such that there exists an auxiliary 
Junction tJj for which one gets, for any positive t, 

^^P{R>x + t^K^)}^^_,_ (1) 

Such a condition imphes that R belongs to the max-domain of attraction of the Gumbel 
distribution; sec Resnick [38], page 26 for more details. De Haan [13] introduced this class 
of distributions as the T-varying class. The function is positive, absolutely continuous 
and satisfies limt^ooV''(0 =0; \iint^ocfp{t)/t = and limt^oo'fpit + xip{t)} /tp{t) = 1 for 
each positive x. It is only unique up to asymptotic equivalence. 

Let us recall Abdous et al. [l]'s result in this rapidly varying context (see Theorem 
l(ii)). 

Theorem 1. Let {X,Y) be a bivariate standardized elliptical random variable with Pear- 
son correlation coefficient p and radial component satisfying (1). Then, for each z G R, 
one has 



lim P{Y <px + z^\- p^yjxtij{x) \X>x)= 4>(z). 

X — >oo 

In the following subsection, a rate of convergence is provided for the approximation 
result stated in Theorem 1, as well as a second order correction. Note that if (X, y) 
has an elliptic distribution with correlation coefficient p, then the couple (X, — F) has 
an elliptic distribution with correlation coefficient —p. Therefore, one can focus on non- 
negative p. From now on, assume that p > 0. As a consequence, one can assume that 
both cc > and y > 0. 



2.1. Main result 

For each distribution function iJ, denote by Pi the survival function H = 1 — H . The 
following assumption will be sufficient to obtain the main result. It is a strengthening of 



Assumption 1. Let H be a rapidly varying distribution function such that 

ii{x + tip{x)} 



H{x) 



for all t>0 and x large enough, where limj;^oo xi^) — Ojip satisfies 

lim = 

K— >oo X 



(2) 



(3) 



and 8 is locally bounded and integrable over [0,oo). 
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The following result is a second order approximation for conditional excess probabilities 
in the elliptical case with rapidly varying radial component. 

Theorem 2. Let {X^Y) be a bivariate elliptical vector with Pearson correlation coeffi- 
cient p Cz [0,1) and radial distribution H that satisfies Hypothesis 1. Then, for all x > 
and z G M, 

p(r <px + z^J\ - p'-^x'i}{x) \ x>x) 

, (4) 



V{Y <px + z^\- p'^^xii{x) + pil}[x) \ X>x) 



= <^{z) + 0\x{x)^ 
All the terms 0{ ) are locally uniform with respect to z. 

Remark 1. This result provides a rate of convergence in the approximation result of 
Abdous et al. [1] Theorem 1, and a second order correction. This correction is useful only 

if x{x)^o{.Ji}{x)/x). 

Remark 2. Theorem 2 and the formula 

f{X<x'] Y<y\X>x) 

= ^{Y <y\X>x)~ P(X >x' \X> x)f{Y <y\X>x') 

yield (with some extra calculations) the asymptotic joint distribution (and the rate of 
convergence) of {X, Y) given that X > x when x is large: for all a; > and z e R, 



''{X<x + tip{x)- Y <px + z^/l- p^y/xi^{x) \ X>x) 



= {l^e-'mz) + o[x{x) + f-^y (6) 
where all the terms 0{ ) are locally uniform with respect to z. 

Remark 3. Hashorva [25] obtained that the conditional limit distribution of Y , given 
that X = X for multivariate elliptical vectors with rapidly varying radial component, is 
also the Gaussian distribution. The equality of these two asymptotic distributions is not 
true in general. In the elliptical context, this is only true if the radial variable is rapidly 
varying. The conditional distribution of Y given that X = a; is of course related to the 
joint distribution of (X, Y) given that X > x via the formula 

P(X < x'; Y <y\X>x)= j T{Y <y\X^ u)Px(du), 
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where Px is the distribution of X . However, the limiting behavior of the integrand is not 
sufficient to obtain the hmit of the integral, so (6) is not a straightforward consequence 
of Hashorva's result. 



2.2. Remarks and examples 

Hypothesis 1 gives a rate of convergence in the conditional excess probability approxima- 
tion. To the best of our knowledge, the literature deals more classically with second order 
conditions providing limits (see, e.g., Beirlant et al. [5], Section 3.3) or with pointwise or 
uniform rates of convergence (cf. De Haan and Stadtmiiller [16] or Beirlant et al. [6]). The 
need here is to have a non-uniform bound that can be used for dominated convergence 
arguments. However, in the examples given below, the non-uniform rates x(a;) that we 
exhibit arc the same as the optimal uniform rates provided by De Haan and Stadtmiiller 
[16]. 

One can, however, check that Hypothesis 1 holds for the usual rapidly varying func- 
tions, in particular, for most of the so-called Von Mises distribution functions, which 
satisfy (see, e.g., Resnick [38], page 40) 



for X greater than some > 0, where d > and ■i/'-ff is positive, absolutely continuous 
and lima;_,oo ?/;^(a;) = 0. Note that under this assumption, -0// = H / H' and is an 
auxiliary function in the sense of (1). In the sequel, auxiliary functions tJjh satisfying (7) 
will be called Von Mises auxiliary functions. 

The following lemma provides sufficient conditions for Hypothesis 1, which could be 
weakened at the price of additional technicalities. 

Lemma 1. Let H be a Von Mises distribution function and assume that 

(i) ipH is ultimately monotone, differentiable and \4''h\ is ultimately decreasing; 

(ii) ifij^H is decreasing, then either VluIx^oo '4' nix) > or there exist positive constants 
ci and C2 such that for all x >0 and u>0, 



■iI)h{x + u) 
Then Hypothesis 1 holds the for ipn '■ 

H{x + til;H{x)} 



< cie"^''. (8) 



e-* 



<xix)eit), 



H{x) 

with x{x) = O [I'll!' fj[x)\) and Q{t) ~ 0{{1 + t)^'^) for an arbitrary k > . 

Remark 4- Assumption (8) holds if V'ff is regularly varying with index 7 < 0. It also 
holds if ipHix) 
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Remark 5. If tpH is regularly varying with index 7 < 1 and "0^ ultimately decreasing, 
then = o{^JiI)h{x) / x) and hence the second order correction is useful. 

Remark 6. This bound corresponds to a worst-case scenario. In many particular 
cases, a much faster rate of convergence can be obtained. For instance, if H{x) = 
e~*, then (2) holds with = 1, whence x(a;) = l/x, but for any positive x and 
t, H{x + tifjuix)} / H{x) = e~*. The rate of convergence is infinite here. If H is the 
Gumbel distribution, then ij^n = 1 and the rate of convergence in (2) is exponential: 
\H{x + t)/H{x) - e-*| < 2e-^e-*. 

Some examples of continuous distributions satisfying Hypothesis 2 are given for the 
purposes of illustration in Section 4, Table 1. The following example illustrates a discrete 
situation. 



Example 1 (Discrete distribution in the domain of attraction of the Gum- 
bel law). Let t/i be a concave increasing function such that lim^j^oo ^(a;) = +00 and 
lim2:_»oo '^''{x) = 0, and define 

iJ#(x)=exp|-£^|, H{x)=H#{[x]\ 

where [x\ is the integer part of x. H is then a discrete distribution function and H 
belongs to the domain of attraction of the Gumbel distribution, but does not satisfy 
condition (7). Nevertheless, following the lines of the proof of Lemma 1, one can check 
that condition (2) holds with 

X(x) ^ 0(^'(N) + l/V-lx)), 9(0 = 0((1 + 

for any arbitrary k> 0. 

One can deduce from Lemma 1 that Theorem 2 holds for ip ~ i^'h ■ The following lemma 
concerns what happens if one uses an asymptotically equivalent auxiliary function 
instead of ipn in Theorem 2. 



Lemma 2. Under the assumption of Lemma 1, let ip be equivalent to TpH at infinity ana 
define ^{x) = \ip{x) — '4'h{x)\/iPh{x)- Then 



H{x + ti'(x)) 



H{x) 



<0{WH{x)\+£,{x)}Q{t)- 



(9) 



Remark 7. A consequence of Lemma 2 is that if an auxiliary function ip is used instead 
of ipH in Theorem 2, then the second order correction is relevant, provided that ^(x) = 
o{\/'Ph{x)/x). This is the case in the examples given on line 1 and line 3 of Table 1 
('Normal' and 'Logis') if one takes ip{x) = 1/x or ip{x) = l/(2a;), respectively, and in the 
example given on line 5 of Table 1 ('Lognor') when using il}{x) — x/log(a;). 
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3. Statistical procedure 

For given large positive x and y, eonsider the problems of estimating 9{x^ y) = P(y < y \ 
X > x) and the conditional quantile function 9{x, •)^. Note, in passing, that in a practical 
situation, x is neither a threshold nor a parameter of the statistical procedure, but a value 
imposed by the practical problem (e.g., a high quantile of the marginal distribution of 
X). The empirical distribution function is useless since there might be no observations 
in the considered range. We suggest estimating these quantities by means of Theorem 2. 

Assume that a sample {Xi , Yi), . . . , (X„, y„) is available, drawn from an elliptical dis- 
tribution with radial component satisfying Hypothesis 2. Note that in this section, we 
do not assume that the distribution is standardized, so a preliminary standardization is 
required. 

3.1. Definition of the estimators 

The estimation of 9{x,y) and the conditional quantile function requires estimates of fix, 
fly, crx, cry, p and ip. Let fix, fJ-Y, o'x, , Pn and ipn denote such estimates. For fixed 
x,y > 0, define 

^„.i(a:, y) = $ f y_^P,^^\ (10^ 



and 



9n,2{x,y) = ^{ I ), (11) 



where x= {x — fix)/cfx and y ~ {y — fiy)/(7y . 

In order to estimate the conditional quantile fimction 0(a;,-)*~, define, for fixed 9 G 

(0,1), 



= /V + o-y{p„i + yjl- pI\J xi^nix)'^ ^{9)}, (12) 

yn,2 = fj-Y + ^yIPuX + Pn^nix) + y/ 1 - pi\J X'f/'n (i)$"^ (6*)} . (13) 

Estimating px , fJ-Y , , o'y and p is a classical topic and the empirical version of each 
quantity can easily be used. Under the assumption of elliptical distributions, however, 
one can observe better stability when the Pearson correlation coefficient is estimated by 
Pn = sin(7T'r„/2), where fn is the empirical Kendall's tau (see, e.g., Hult and Lindskog 
[30] for more details). This estimator is y^-consistent and asymptotically normal. 

Consider now the problem of the estimation of ip. Since auxiliary functions are defined 
up to asymptotic equivalence, a particular representantive must be a priori chosen in 
order to define the estimator. We assume that an admissible auxiliary function is 
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i>{x) = ^x'-^ (14) 

for some constants c > and /3 > 0. Under this assumption, estimation of ip reduces to 
estimating c and f3. 

An extensive body of literature exists on estimators of /? (see, e.g., Bcirlant et al. [4], 
Gardes and Girard [23], Dicrckx et al. [17], among others). The method chosen here is the 
one proposed in Bcirlant et al. [7]. Let A:„ be a user-chosen threshold and \ < j <n, 
be the order statistics of the sample i?i , . . . , i?„ . The estimator of (3 is obtained as the 
slope of the WeibuU quantile plot at the point (loglog(n/fc„), log(i?„_fc^_„)): 



p K ^ T.'iZi log ^ogjn/i) - log logjn/kn) ^^^^ 

kn^ Y!IZi l0g(i?„_i+i,„) - l0g(i?„_fc„,„) 



An estimator of c is then naturally given by 



1 ^ log(n/i) 



■ 

Actually, in our context, the radial component is not observed. We estimate the i?i's by 

where Xi = [Xi — fix)/(^x and Yi = (Yi — fi,Y)/^Y, and we plug these values into (15) 
and (16). We then define i^nix) = a;^~^"/(c„/3„). 

3.2. Discussion of the estimation error versus approximation 
error 

Let 

_ y- PnX , _ y - PnX - Pnlpnjx) 
^n,l — I ; ^n,2 — I 1 

V^~PI\ ^'^n (x) V^-Ply ^'^n (x) 

y — px y — px ^ px^^^ /{c(3) 



y^{r^y/x^-p/{cpy ^/{r^^X^~P/{c(3) 

Then, for i = 1,2, 

dn4^, y) - 0{x, y) = $(z„,,) - $(z,) + $(zO - e{x, y). 

This shows that the estimators defined in (10) and (11) have two sources of error: the 
first one, ^{zn^ — '^{zi)^ comes from the estimation of p, p, a and -0, and the second 
one, ^(zi) — 9{x,y), from the asymptotic nature of the approximations (4) and (5). 
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The order of magnitude of the estimation error can be measured by the rate of con- 
vergence of the estimators. In order to obtain a rate of convergence for the estimators 
/3„ and c„, we assume that H is a, Von Mises distribution function with 

^H{x) = \x'-P{l + t{x)}, (17) 

where i is a regularly varying function^ with index ?7/3 for some jy < 0. This implies that 
H{x) ~ exp{— ca;^[l + s(x)]}, where s is also regularly varying with index -qP. Under this 
assumption, the function ip defined in (14) is an admissible auxiliary function and Girard 

^ 1/2 

[24] has shown that fin is /c„ -consistent, for any sequence fc„ such that 
fc„^oo, A:y2iog-i(n/fc„)^0, kl'H{\og{n/kn)) ^ 
where b is regularly varying with index 77; see Girard [24], Theorem 2 for details. Similarly, 

1 /2 

it can be shown that fc„ (c„ — c) = Op{\) under the same assumptions on the sequence 

1/2 

kn- Thus, for any x, tpn{x) is a kn -consistent estimator of tp{x). Besides, p„, fin and an 
are -v/n-consistent, so the estimation error ^{zn^i) — ^{zi) is of order kn in probability. 

If if>H satisfies (17), then Hypothesis 1 holds, and Theorem 2 and Lemma 2 provide a 
bound for the deterministic approximation error $(2^) — 9{x,y). Some easy computation 
shows that the second order correction is useful only if 77 < —1/2. 



3.3. Discussion of an alternative method 

The method described previously makes use of the asymptotic approximations of Theo- 
rem 2. It could be thought that a direct method, making use of (19) and an estimator of 
H, would yield a better estimate of 9{x, y) since it would avoid this approximation step. 
Recall, however, that we specifically need to estimate the tail of the radial distribution 
so that a nonparametric estimator of H cannot be considered. The traditional solution 
given by extreme value theory consists of fitting a parametric model for the tail. This 
will always induce an approximation error. 

Nevertheless, there exists a situation in which the approximation error can be canceled: 
if the radial component is exactly WeibuU distributed, that is, H{x) = exp{—cx^}, then 
tpH satisfies (14), so for any a;,u > 0, (7) implies that 



H{xu) 




Therefore, in this specific case, a consistent estimator of 0{x,y), say On.3{x,y), can be 
introduced which does not make use of any asymptotic expansion as in Theorem 2. Via 

function / is regularly varying at infinity with index a if for all t>0, limx^oo f{t^)/f{^) = t"- 
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Table 1. Bivariate elliptical distributions used for the simulations, together with elliptical gen- 
erator g, Von Mises auxiliary function ipH (or equivalent), the function x defined in Hypothesis 1 
and the values of the parameters used (in addition to p G {0.5,0.9}) 

Bivariate law Generator g{u) ipH{x) x(^) Parameters 



Normal e""'^ 

Kotz u''/'-'e-"''^' x^-^/P 0{x-^) /?G{1,4} 

T • * e-" 1 + e-"' ^, _2, 

ff*(«)* ^T-^^T^ Of-^) /3 = 3/2 

^ ^ l + /?loga:: Vlogs/ ' 

Lognor" ie-('-^")/« _^ + of^) of ^ 

M logs Vlog'^a;/ Vlog 



Modified Kotz 



Student 



log"^ xJ V log s / 

fl + ^j - - i^G{3,20} 

* "Logis" and "Lognor" refer to the elliptical distributions with logistic and lognormal generator, respec- 
tively. 

**g*(«) = {(3/8) logn + l/2}n-i/4 exp{-(l/2)«3/4 iog„}. 



(19), we get, explicitly. 



'^'^ k{x, X, cos(u)) All + r^*^'^^*") K{x, y, sm{u + Uo)) du 



> arctan(ro) \ ' ' \ / / J — jj^ 

2 J^^'^ K{x,x,cos{u))du 



V3(a;,y)^l- "— ^-"^ _ , (18) 



where 



K{x, y, v) = exp 



ds 



X lpn{s) 



io = [y/x - /0„)/ Vl -/5^, C^o = arctan(/5„/ ^1-/5^), 
£ = (a; - /ix)/o-x, y ^ {y - fiy) / cry ■ 

We included this estimator in the simulation study as a benchmark when looking at 
elliptical Kotz-distributed observations (see Table 1). 



4. Simulation study 

To assess the performance of the proposed estimators, we simulated several families of 
bivariate standard elliptical distributions. Recall that a standardized bivariate elliptical 
density function can be written as f{x,y) = Cg{{x^ — 2pxy + ?/^)/(l — p^)}, where g 
is called the generator, p is the Pearson correlation coefficient and C is a normalizing 
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constant. The density of the radial component R is given by H'{r) = Krg{r^), where K 
is a normahzing constant (see, e.g., Fang et al. [21]). 

The distributions used are presented in Table 1. The Pearson correlation coefficient 
will be either p = 0.5 or p = 0.9. Three of them (Normal, Kotz and Logis) are Von 
Mises distributions which satisfy both Hypothesis 1 and (17); in addition, the Von Mises 
auxiliary function of the Kotz distribution satisfies (14). The Lognor and the modified 
Kotz distributions satisfy Hypothesis 1 but not (14) and, finally, the bivariate Student 
distribution has a regularly varying radial component, so it does not satisfy any of the 
assumptions. These three distributions arc used to explore the robustness of the proposed 
estimation method. 

In each case, 200 samples of size 500 were simulated. Several values of x were cho- 
sen, corresponding to the (1 — p)-quantile of the marginal distribution of X , with 
p = 10"'^, p = 10~^ and p = 10~^. For each value of x, we computed (by numeri- 
cal integration) the theoretical values of y corresponding to the probability 0{x,y) = 
0.05,0.1,0.2, . . . ,0.8,0.9,0.95. Wc then estimated 0{x,y) via the three proposed methods 
(cf. Section 3). For the estimation of the auxiliary function ip, the threshold chosen cor- 
responds to the highest 10% of the estimated Ri's. It must be noted that this choice is 
independent of x. For each fixed x, we also estimated the conditional quantile function 
9{x,-)^ by both methods (12) and (13). We did not compute the estimated quantile 
function for Method 3 since it would involve the numerical inversion of the integrals 
which appear in (19). This is one advantage of Methods 1 and 2 over Method 3. 

Some general features can be observed which conform to theoretical expectations, (i) 
First, in the given range of x and y, there were hardly any observations, so the empirical 
conditional distribution function is useless, (ii) For a given probability 6, the variability 
of the estimators slightly increases with x for all underlying distributions. For a given x, 
the variability of the estimators is greater for medium values of 9. (iii) The results for the 
Student distribution are as expected: if the degree of freedom v is large, the estimation 
shows a high variability but moderate bias, while if v is small, then the estimation is 
clearly inconsistent, (iv) The results for the Logis and modified Kotz distributions are 
similar to those for Gaussian distribution, (v) As described in Section 3.3, Method 3 is 
markedly better for the Kotz distribution only. 

Hence, wc have chosen to report only the results for the largest value of x (correspond- 
ing to the 10^^-quantile of the marginal distribution of X) and the Normal, Kotz (with 
parameter P = I and /3 = 4) and Lognor distributions, for p ~ 0.5 and p ~ 0.9. Figures 1- 
4 illustrate the behavior of the estimators of the probability 9: median, 2.5% and 97.5% 
quantilcs of the estimation error 9„^i{x,y) — 9{x,y) {i — 1,2,3) are shown as a function 
of the estimated probability. Figure 5 shows the estimated conditional quantile functions 
yn,iix, y) {i = 1, 2) and the theoretical conditional quantile function y = 9{x, •)^ for only 
three distributions and p ~ 0.9 because the results are much more stable as the correlation 
or the distribution vary. Median, 2.5% and 97.5% quantilcs of the estimated conditional 
quantile function y„^i(x,y) (i = 1,2) are given as a function of the probability. 

From these simulation results, one can see that the estimator of 6 by Method 1 presents 
a systematic positive bias which, of course, induces an underestimation of the conditional 
quantile function. As expected. Method 2 corrects this systematic bias; the correction is 
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Figure 1. Median, 2.5% and 97.5% quantiles of tlie estimation error 6n,i{x,y) — 6{x,y) 
(i = 1, 2, 3) as a function of tlie estimated probability. Gaussian distribution. First row: p — 0.9; 
second row: p = 0.5. 



better when p is large. This is also true for the Lognormal generator, though to a lesser 
extent. 

As already mentioned, the Lognor and modified Kotz distributions do not satisfy the 
assumption (17). In both cases the radial component belongs to an extended WeibuU-type 
family, with auxiliary function ijj of the form 

V'(x) = cx^-^{logx)-^{l + o(l)} 

with c > 0, /? > and S > if /3 = 0. The modified Kotz distribution corresponds to 
(3 = 3/2 and 6 = 1 and Lognor corresponds to /? = and 6 = 1. The simulation results are 
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Figure 2. Median, 2.5% and 97.5% quantiles of tlie estimation error 6n,i(x,y) — 6{x,y) 
(i — 1, 2, 3) as a function of tlie estimated probability. Kotz distribution, /3 = 1. First row: p = 0.9; 
second row: p = 0.5. 



much better for the modified Kotz than for the Lognor distribution. This tends to prove 
that the method is not severely affected by the logarithmic factor, as long as /3 > 0. 

5. Financial application 

In this section, the practical usefulness of our estimation procedure is illustrated in the 
context of financial contagion, for which an estimation of the conditional excess proba- 
bility is needed. Data used by Levy and Duchin [35] are here revisited. More precisely, 
we consider the series of monthly returns for the 3M stock and the Dow Jones Industrial 
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Figure 3. Median, 2.5% and 97.5% quantiles of tlie estimation error 9n,i{x,y) — 6{x,y) 
(i = 1, 2, 3) as a function of tlie estimated probability. Kotz distribution, /3 = 4. First row: p = 0.9; 
second row: p = 0.5. 



Average for the period ranging from January 1970 to January 2008. In the sequel, we 
arbitrarily investigate the conditional behavior of the 3M stock monthly returns, given 
some extreme values of the Dow Jones Industrial Average. 

According to Levy and Duchin [35], these two series can be marginally fitted by a 
logistic distribution. Indeed, a Kolmogorov-Smirnov goodncss-of-fit test of the logistic 
distribution gave us a P-value of 0.48 for the 3M returns and 0.49 for the Dow Jones 
Industrial Average returns. These P-values were obtained via a Monte Carlo simulation, 
following the procedure outlined by Stephens [40]. Moreover, the test of elliptical sym- 
metry of Huffer and Park [29] was used to show that the data fit the bivariate elliptical 
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Figure 4. Median, 2.5% and 97.5% quantiles of tlie estimation error 9n,i{x,y) — 6{x,y) 
(i — 1,2,3) as a function of the estimated probability. Lognor distribution. First row: p = 0.9; 
second row: p = 0.5. 



model (P-value = 0.61). Finally, we checked that both marginal upper tails exhibit rapid 
variation; for this, we performed a generalized Pareto distribution fit to the 15%-largest 
values and checked via a test based on the profile likelihood 95%-confidence interval that 
the shape parameter could be considered equal to at level 5% (see, e.g.. Coles [12] for 
details on these classical procedures). Note that the estimated generalized Pareto tail 
agreed completely with the tail of the fitted logistic distribution in both cases. 

Consequently, the estimation procedures presented in the previous sections can be 
applied to these data. As an illustration, in Figure 6, we depict the three estimates of y i-^ 
1 — 9{x, y) = f{Y > y\X > x) for different values of x corresponding to the 0.975, 0.99, 
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Figure 5. Median (solid line), 2.5% and 97.5% quantiles (dashed lines) of the estimated con- 
ditional quantile function yn,i {i = 1,2) defined in (12) and (13) and theoretical conditional 
quantile function y (dotted line) as a function of the probability 9. First row: Normal distri- 
bution; second row: Kotz distribution, /3 = 4; third row: Lognor distribution. For each of these, 
p = 0.9. 



0.999 and 0.9999 quantiles of the fitted logistic distribution, together with the estimated 
marginal survival function F{Y > y). This last probability was estimated via the logistic 
distribution fitted to the K^'s. It is clearly evident from these graphics that 1 — ^„_2 and 
1 ~ ^n,3 provide very similar estimates, whereas 1 — 9„^i gives uniformly smaller values. 
All these estimates are uniformly greater than the marginal survival function of Y. This 
allows us to conclude that the data exhibit contagion from the Dow Jones Industrial 
Average to the 3M stock. 
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Figure 6. Estimates y i— > 1 — 9n,i{x,y) [i — 1,2,3) of tlie conditional excess distribution of the 
3M stock monthly return, given that the Dow Jones Industrial Average monthly return exceeds 
four extreme values. The subplots (a), (b), (c) and (d) are for the values x corresponding to 
the 0.975, 0.99, 0.999 and 0.9999 quantiles of the fitted logistic distribution, respectively. The 
solid line is for 0„^i, the dashed line for 8„.2 and the dotted line for ^„,3. The estimate of the 
marginal survival function P{Y > y) is shown as a dotted-dashed line. 



6. Concluding remarks 

In this paper, we restricted our attention to elliptically distributed random pairs {X, Y) 
having a rapidly varying radial component. Three methods have been proposed to esti- 
mate the conditional excess probability 6{x,y) for large x. Under this specific assumption 
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of elliptical distributions, Methods 2 and 3 revealed comparable results and outperformed 
Method 1. 

Methods 1 and 2 make use of an asymptotic approximation of the conditional excess 
distribution function by the Gaussian distribution function. As shown by Balkema and 
Embrechts [3], this approximation remains valid outside the family of elliptical distri- 
butions, under geometric assumptions on the level curves of the joint density function 
of {X, Y) . This suggests that these methods may be useful outside the elliptical family. 
This is an ongoing research project. 



Appendix 



Proof of Theorem 2. Define Uq = arctan(/9/-\/l — p'^). For a; > and y £ (0,a;), we 
have 



fix >x,Y >y) 

X 

cost/ " pcos[/+ -P^sinC/ ' """""2 



R>-^V ^ ; -Uo<U<^ 



Let to ~ {y/x — p) / V 1 ^ P^- Then — [/q < arctan(io), and x/cosu > y/{pcosu + 
\J\ — p^ sinu) if and only if w > arctan(to)- Hence, 

Y{X>x,Y >y) 



-/2 - d^Z ta„(to) r y ^ 



' arctan(iQ) 

p(r > y I X > x) 



cosM J 27t J -Va \ sin(u + [/q) J 27t' 



dw 



2 /J'^^i/(x/cos(w))dw 



(19) 



If to > 0, that is, ?/ — > 0, the changes of variables v =1/ cos(ti) and v = 1/ sin(u + Uq) 
yield 

P(y > y I X > .t) = :^L±ii (20) 



with 



H(vx) dv 



■1 i?(x) t,V^;2^' 
H{vy) dv 
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H{vx) Av 



h = 2 



1 H{x) t;v^;2— T' 



W2 = xwi/y. 



If io < 0, then 



nY>y\X>x) = ^' (21) 
Denote wq = x{wi — l)/^p{x). In Ii and I3, the change of variable v = 1 + -^^^ t yields 

h - 



i'jx) H{x+tipix)) dt 

J Wo 



{l+tiiix)/x))^l + it/2)ii,ix)/x)V2't 



J ^2 /^(^) rH{x+t^{x)) dt 



(l + t(V'(a;)/x))^l + (t/2)(^(:r)/x)y2t 
In the change of variable vy = x -\- t'ip{x) yields 

, rnix + t^ix)) iy/x)dt 

h = 



Hix) (1 + t{4,{x)/x))^l - {y/xY + 2t{i,{x)lx) + {^^^ {x) / x'^)f^ ' 

Let Ji = ^Jx/'4>{x)Ii, i = 1, 3, and J2 = {x/ip{x))l2- We start with /i and I3. We will use 
the bound, valid for all B,C > 0, 

0<1-- ^ , <B/2 + C, (22) 



which follows from straightforward algebra and the concavity of the function x ^ \Jl + x. 
Applying this bound with B = and C ~ ^^-^t yields 

- ~ (1 +t(V^(x)/x)) VI + (i/2)(V'(^TO " 4^ 

We thus have 



I Ji - V2n<t>{V2uJo)\ + \J3- V2n\ 
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-^=$(^/2«^) + 0(^x(x) + 
Now consider J2. Applying the bound (22) with C = tip{x)/x and 
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(23) 



B=l2t 



and making use of Hypothesis 1, we obtain 
[y / x)\/2jiip{^/2w^) 



J 2 



< 



y/x 



X{x) / Q{t)dt + 



y/x ip{x) 
1 — {y/xY X 



ip{x) 



Choose y = px + ^ \/ xtJj {x)z for some fixed z G M. Then, for large enough x, we 
have <y < X and 

y/x p 



-0(VV<^) 



Thus, 



k 



p 



%Ij{x) 



V'(x) 



(24) 



For z > and large enough a;, plugging (23) and (24) into (20) yields 

P 



0(x,2/) = $(V2^) 



i}{x) 



(/9(x/2^)+ O x{x) + 



iix) 



^Jl-P^ V X 

For z < and large enough x', plugging (23) and (24) into (21) yields 



0(a:,2/) = l'(V2^) 



?/>(x) 



(^(V2^)+ O x(a;) + 



i^{x) 



Now note that ujq = / 2 + 0{il}{x) / x) , hence \/2wq = \z\ +0{ip{x)/x). Thus, in both the 
cases z >0 and z < 0, (4) holds. Let z = z' + py/ ip{x) /x/ y^l — p^. A Taylor expansion 
of $ and ip around z' yields (5). □ 
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Proof of Lemma 1. 



H{x + tiP{x)} 

- — e = exp 



H{x) 



cxp 



■d.s 



/o V'{a; + sV(a;)} 

Applying the inequality |c~° — e~^| < |a — 6|e^°^^ valid for all a, 6 > yields 
H{x + til}{x)} 



< 



ds 


exp 


-tA [ 






Jo 



Let 



and 



Hix) 

* Vi^; + 3^/^(0;)} — tp{^) 

* l^pjx + S'ijjix)} - tlj{x)\ 
ip{x + stjj{x)} 

i;{x) 



ij{x) 



ip{x + si'{x)} 



ds 



ds = I{x, t) 



exp 



-i A 



Tp{x + si}}{x)} 



ds 



-■E{x,t). 



Case tp increasing. If ip is non-decreasing, then t/"' > and -0' is decreasing. Thus, for 
any i5 > and large enough x, ip'ix) < 6 and 



ip{x + sip{x)} Jq tJj{x) + s^' {x)'ip{x) 



ip{x) 



■ds 



ds 



> 



* ds 1 

■7\og{l + St). 



1 + stJj'{x) - Jf^ 1 + sS S 

This implies that E{x,t) < (1 + 6t)^^/^ for large enough x. Since -0 is increasing and ip' 
is decreasing, we also have 



— ds < lb ix) — . 



Thus, for any 5 >0, l{x^)E(^x,t) =0{\i\)' (x)\t^(\ ^ 8t)-^l^). 
Case decreasing. If ip is monotone non-increasing, then 



ip{x) 



Tp{x + S1p{x)} 



ds>t 
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and E{x,t) <e~*. Also, since \iJj'\ is decreasing, 

I{x,t)<W{x)\ 
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sip{x) 



tp{x + sip{x)} 



ds. 



If ip has a positive limit at infinity, then I{x, t) = 0{ip' {x)t'^). Otherwise, lim^j-Kxa ''Pi^) = 
and (8) holds. This yields, for large enough x, 



I{x,t)<\tp'{x)\ 



sip{x) 



ipix + stJj{x)} 



ds 



<ci\^'ix)\ / se'=='''^(^Ms<ci|V''(a;)| 



'/^ds. 



Thus, I{x,t)^Oi\^'{x)\t^e*/^) and I{x,t)E{x,t) = 0{\'ilj'{x)\t^e-*/^). This concludes the 



proof. 

Proof of Lemma 2. Write 
H{x + tip{x)} 



□ 



H{x) 



Hix + tipix)} 



< 



H{x + tilj{x)} 



H{x + tipHix)} 

H{x + tip{x)} 
H{x + tipH{x)} 

tl])(x)/lpH(x) 



H{x + tipH{x)} 



exp 



H{x + tipHix)} 
If ^jh is increasing, then 

ftlf}{x)/lf!Hix) 



H{x) 
^h{x) 



lj}{x + si!H{x)} 



ds 



ipnix) 



exp 



■iPh{x + sipnix)} 



ds 



ds 



< ti{x) 

< e*«(") 



^pH{x + s^pHix)} 
Since £,{x) 0, gathering (27) and (28) yields, for large enough x, 

Hix + tipix)} 



H{x + tipHix)} 



1 



e-*<t^(x)e-*/2. 



(25) 
(26) 

(27) 
(28) 

(29) 



If TpH is decreasing and lmix^oo4'Hix) > 0, then the ratio '4>h{x)/'iPh{x + sTp{x)} is 
bounded above and away from 0, so 



ttl>{x)/4:H{x) 



iPh{x) 



ipix + sipnix)} 



ds 



<ctax), 



and (29) still holds. 

If ipH is decreasing and lim^j^cx) i^nix) ~ 0, then applying (8) gives that the left-hand 
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side of the previous equation is bounded by Ci^(x) eyi]i{2c2ip{x)i\ . Thus, for large enough 
X, (29) still holds. 

This provides a bound for the right-hand side of (25). The term (26) is bounded by 
Lemma 1. □ 
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